Pairing states of a polarized Fermi gas trapped in a one-dimensional optical lattice 
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We study the properties of a one-dimensional (ID) gas of fermions trapped in a lattice by means of the density 
matrix renormalization group method, focusing on the case of unequal spin populations, and strong attractive 
interaction. In the low density regime, the system phase-separates into a well defined superconducting core 
and a fully polarized metallic cloud surrounding it. We argue that the superconducting phase corresponds to 
a ID analog of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, with a quasi-condensate of tightly bound 
bosonic pairs with a finite center-of-mass momentum that scales linearly with the magnetization. In the large 
density limit, the system allows for four phases: in the core, we either find a Fock state of localized pairs or 
a metallic shell with free spin-down fermions moving in a fully filled background of spin-up fermions. As the 
magnetization increases, the Fock state disappears to give room for a metallic phase, with a partially polarized 
superconducting FFLO shell and a fully polarized metallic cloud surrounding the core. 
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Ultracold atoms in optical lattices can be used to study 
models of strongly correlated fermions in clean and controlled 
experimental conditions. In particular, cold gases provide an 
optimal playground to study the crossover between a Bardeen- 
Cooper-Schrieffer (BCS) superfluid, with extended Cooper 
pairs, and a Bose-Einstein condensate (BEC), composed of 
molecules of tightly bound pairs. 1 As cold atom gases can 
also be realized in optical lattices, dimensional crossover ef- 
fects can be accessed. 2 In particular, ID optical lattices can 
be prepared by strongly ramping up the amplitudes of two 
out of three counterpropagating light waves. 2 This, combined 
with the possibility of tuning the interactions, allows for the 
realization of the fermionic ID Hubbard model, with, in 
experiments, two hyperfine states interacting via an onsite 
potential. 3 In the case of a spin imbalanced fermion mix- 
ture, magnetized superconducting states are expected, such 
as, e.g., a superfluid-normal mixture 4 or the one predicted by 
Fulde and Ferrell 5 and Larkin and Ovchinnikov 6 (FFLO state) 
over four decades ago. The FFLO state 7 is characterized by 
pairing across a spin-split Fermi surface, with the resulting 
Cooper pairs having a finite center-of-mass momentum, pro- 
portional to the spin polarization, and consequently, an oscil- 
latory phase in the superconducting correlation function. Its 
experimental observation has eluded condensed matter physi- 
cist until very recently, when it was detected in heavy-fermion 
systems. 8 Its realization in cold atom systems has acquired 
particular relevance for the field of high-T c superconductiv- 
ity. Moreo and Scalapinc 9 have recently pointed out that, by 
exploiting a particle-hole transformation, 10 its presence in the 
2D attractive Hubbard model may imply the existence of a 
striped phase in its repulsive counterpart. 

Experiments on 3D trap a 11 ' 12 suggest the existence of 
phase-separated shells, with a superfluid core, and a par- 
tially polarized normal cloud, and a future characterization of 
the superconducting phases may establish whether the FFLO 
phase is present or not. Theoretical work indicates that the 
FFLO phase is stable in a narrow window around the (unpo- 



larized) Fermi surface. 7 However, the instability against an 
FFLO state may be enhanced in low dimensions, rendering 
this state more robust. 

The uniform ID polarized Fermi gas has been stud- 
ied by means of bosonization and renormalization group 
techniques, 13 providing evidence for the existence of an FFLO 
state. The case of fermions confined to a ID trap has 
been analyzed using modified versions of the Gaudin-Yang 
Hamiltonianr 4 ^ a minimal integrable model of fermions in 
the continuum interacting via a contact potential. A two- 
shell structure has been predictedr^ 7 - with a partially polar- 
ized phase of the FFLO type in the center of the trap, and 
either fully paired or fully polarized wings, depending on the 
strength of the magnetic field or, equivalently, the total mag- 
netization. 

In this work we investigate the FFLO state in a Hubbard 
chain, thus accouting for the optical lattice: 

L-l L 

H = - t (c\ a c i+lrJ + h.c^j + U ^2 n^riii 

i—\,a i—1 
L 

+ vY / ^~L/2) 2 n l , (1) 

2 = 1 

where ct creates a fermion with spin a =|, j at site I, 
nio = electa, ni — n/-\ + ngi is the local density, t is the 
hopping parameter, and U is the onsite interaction energy, 
which in this work is negative. We define x = ia, where a 
is the lattice spacing, set to unity. We add a harmonic con- 
fining potential parameterized by a constant V . The Hub- 
bard model with V — has been extensively studied, and its 
properties are well documented in the literature. 18 The low- 
energy properties of the Hamiltonian (Q]) with U < are 
those of a Luther-Emery liquid. 19 At small attractive interac- 
tions, fermions form Cooper-pair-like bound states with a spin 
gap, reminiscent of the superconducting gap in conventional 
BSC superconductors. In the case of strong interactions, the 
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FIG. 1: (color online) Results for the ID attractive Hubbard model 
(U = — St, L = 80), with N = 40 fermions confined by a parabolic 
potential with strength V = 0.002t, for different magnetizations: 
(a) Density profile, (b) magnetization profile, (c) charge fluctuations, 
and (d) spin fluctuations. 



pairs become tightly bound with their extension of the order 
of the lattice spacing only, effectively behaving as hard-core 
bosons, while the spin degree of freedom moves to high ener- 
gies. Superconducting correlations then decay algebraically. 
The parabolic trap adds an extra ingredient that leads to the 
emergence of non-uniform states, and the inclusion of a lat- 
tice may produce a richer phase diagram. It has been shown, 
for instance, that this allows for the possibility of engineer- 
ing states of hard-core bosons such as pure Fock states that 
cannot be found in the continuum. 20 At the same time, the 
physics of the system without a lattice can be recovered in 
the low density limit. While recent studies have addressed the 
unpolarized situation 21 , this article's chief case is the polar- 
ized one. Note that the corresponding situation with repulsive 
interactions has been studied in, e.g., Ref.l22l 

We use the density matrix renormalization group (DMRG) 
method 23 to obtain the ground state properties of this model in 
finite systems, for different numbers of particles N ~ + 
N±, and values of the total magnetization S z = (iVf — -/VjJ/2. 
In our calculations we choose L = 80, U = —8t, and V — 
0.002i, unless otherwise stated. In the following we describe 
two typical situations in the low and large density regimes, 
which illustrate the main features of the problem. 

We start by looking at the low density case with N = 40 
particles. In FigQ]we show the local density (rii) and the spin 
projection (Sf) along with the fluctuations of both quantities. 
We define the charge fluctuations as ((Srii) 2 ) = (nf ) — (rii) 2 , 
and similar for the spin fluctuations ((SSf ) 2 ). 

The system exhibits a nearly half-filled density profile in 
the center with a sharp edge, and oscillations accompanied 
by large fluctuations. As the magnetization grows, the core of 
the system becomes partially polarized, displaying charge and 
spin oscillations, surrounded by fully polarized clouds with 
spin f-fermions only. We can clearly identify two well de- 
fined phases: the fully polarized metallic wings, surrounding 
a polarized state in the core that we wish to characterize in 
more detail in the following. 
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FIG. 2: (color online) (a) Pair momentum distribution function for 
N = 40 fermions with iV pair = J2 k n T" ■ The dotted lines snow 
results for the sum in Eq. (O restricted to the center of the trap, 
(b) Spatial decay of pair correlations (squares). Dashed line: Fit to 
n im' r 0C cos (k \ l — m\) / \l — rn\ a (compare Ref. 13). Data shown 
for L = 100, N = 80, V = 0.0005£, S z = 8, corresponding to the 
same effective density p c ff = NyV '. (c) Momentum fcmax, at which 
the distribution shown in (a) is peaked, vs. magnetization of the core 
S* s /L cB (L cff » 40, see text). 

For this purpose, it is illustrating to analyze the natural pair 
excitations of the system in terms of a bosonic description 
with creation and annihilation operators b\ = cL cj. and 6j = 
Cj^Cif. In the low-density limit, indeed, the pairs behave in 
good approximation like bosons since [b\, bi] = 1 — 2rii m 1, 
but we will use it generically, as it provides an intuitive pic- 
ture. We further define the pair one-particle density matrix 
(OPDM) as Pij = (b\bj). 

In Fig. |2] we show the pair momentum distribution function 
(MDF) 

^pair = Qj L j ^2 exp[ik(l ~ m)] p im . (2) 

l in 

In the unpolarized case, we see a sharp peak centered at mo- 
mentum k = 0, as observed for quasi-condensates of hard- 
core bosons. 20 As the magnetization increases, the distribu- 
tion exhibits two maxima, centered at momenta ±& max . Note 
that by restricting the summation over l,m to the core part 
of the systems, the features in n^ n become more prominent 
(dotted lines in FigO. From Fig. [2] (c), we see that mo- 
mentum fc max is proportional to the magnetization in the core 
S e s, that is obtained by integrating (Sf) over the effective 
size of the region occupied by the FFLO state, from FigQ](a). 
Our results reproduce precisely the behavior expected for the 
FFLO state, in which the pairs possess a finite center-of-mass 
momentum k = kp-\ — kpi (where kp a is the Fermi vec- 
tor of the spin-er fermions), which is predicted to grow as 
k = 7r(n T - ni)/L = 2ttS z /L. 7 In Fig]2](b) we show the 
spatial decay of the pair correlations, consistent with a power- 
law decay of the form nf™ 1 oc cos(/c max |Z — — m\ a , as 
predicted in Ref. [l3i 

The bosonic quasi-condensate can be studied by means 
of Penrose and Onsager's description of the superfluid or- 
der parameter^ The natural orbitals (NO) ip a of the system 
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FIG. 3: (color online) Natural orbital of the pair OPDM for N = 
40 confined fermions and different magnetizations. The inset shows 
the OPDM's eigenvalues A Q 
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FIG. 4: (color online) Same parameters as in FigQ] but N = 80 
fermions. (a) Density profile, (b) magnetization profile, (c) charge 
fluctuations, and (d) spin fluctuations. 



will simply be the single particle eigenstates - in the bosonic 
sense - of the pair OPDM, and the corresponding eigenval- 
ues X a represent their occupations. The NO with the largest 
eigenvalue, ip , is the single-particle state in which quasi- 
condensation takes place. The lowest NO and the eigenvalues 
of the OPDM are presented in the main panel and the inset 
of Fig. |2 respectively. The distribution of eigenvalues qual- 
itatively resembles the results for a trapped gas of hard-core 
bosons. 20 In the unpolarized case we observe a macroscopic 
occupation of the lowest eigenstate, corresponding to the or- 
der parameter. The effect of increasing the magnetization is to 
decrease the occupation of the bosonic condensate, rendering 
the profile of A Q vs a less pronounced. The lowest NO, rep- 
resenting the bosonic order parameter, resembles the density 
profile in the unpolarized case, see Fig. [Tf a )- As the mag- 
netization increases, it develops sharp oscillations, in agree- 
ment with the results for the MDF. As observed in Ref. |s£ we 
see that the unpaired fermions accumulate in the nodes of the 
order parameter, effectively forming magnetic domain walls. 
The single-particle wave function is well confined to the core 
of the system, an indication of the phase separation between 
the FFLO state in the center and the fully polarized metallic 
wings. This observation, i.e., the confinement of the NO to a 
window of length L c g ~ L/2 further corroborates the use of 
the effective quantities L c s and S* s in the inset of Fig.|2jb). 

We now turn our focus to the dense limit, by looking at a 
system with N = 80 particles. The results for the local den- 
sity and magnetization are shown in Fig. [4] along with their 
fluctuations. In the unpolarized case S z = 0, a plateau at den- 
sity (rij ) — 2 appears in the center, surrounded by an extended 
region with oscillations around (rii) = 1. The state in the cen- 
ter is a pure Fock state of localized pairs, without coherence, 
and decoupled from the rest of the system. As a consequence, 
the charge and spin fluctuations in the center of the trap are 
totally suppressed, which are signatures of a band insulator 
with localized bound pairs, and both finite charge and spin 
gaps. The fluctuations increase abruptly at the boundary with 
the partially-filled, unpolarized region. As the total magneti- 
zation increases, pairs are broken and the polarized fermions 
move toward the edges, effectively suppressing the gaps and 



the band insulating behavior in the center of the trap. The 
Fock state survives small values of S z but disappears rapidly, 
and the local magnetization becomes almost uniform, roughly 
equal to the average total magnetization, (S?) « S z /N. The 
fluctuations are most prominent in the intermediate region be- 
tween the wings and the core, and after reaching a maximum 
there, they decrease again to values similar to those observed 
in the center. 

More details are revealed in Figs. [2 a) and (b), where we 
plot the densities (rii a ) for S z = 2 and 8, respectively. Ob- 
viously, the wings are fully polarized, and the f-spins fill the 
core of the trap with a plateau of density (n^) = 1. The Fock 
state survives at S z = 2 and we find a distribution (riii) that 
resembles the one of trapped spinless fermions. 20 Fig. [3c) 
shows the MDF calculated in shells centered in the middle of 
the trap, by restricting the sum in Eq. (O- As we cross the 
different phases surrounding the core, the distribution evolves 
from the bulk result with two peaks to a featureless uniform 
profile describing localized pairs. This shows that in the cen- 
tral region, only the .[-fermions exhibit any dynamics. 

These results establish the existence of three distinct 
phases: first, there is a metallic, fully polarized phase at 
the wings, behaving like a partially-filled band of spinless 
fermions and second, a partially polarized phase in the center 
that is a metal of freely moving j-fermions in a uniform back- 
ground of f -fermions, with an effective site energy of — U. An 
intermediate, third shell separates the two, with features that 
resemble those observed in the FFLO state in the low den- 
sity regime. Notice that a fourth phase, the Fock state in the 
center, survives at small values of the magnetization, before 
giving way to the metallic core. 

In order to shed light on the nature of this intermediate 
phase, we analyze the pair OPDM's spectrum and NOs. In 
Fig.|5jd), we plot the first three NOs that do not correspond to 
localized bound states. These NOs are precisely situated in the 
region where the intermediate phase is found. We have also 
calculated the MDF by restricting the summation in Eq. (0 to 
only this region of the lattice (and the reflected sector on the 
opposite side). This restricted MDF exhibits finite center-of- 
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FIG. 5: (color online) (a),(b) Particle density for (n; CT ), a =\, J., 
for TV = 80 fermions, and S z = 2, 8, respectively. The double 
occupancy (di) = (n^nn) is also included, (c) Pair momentum 
distribution function measured in concentric shells of different radii, 
for S z = 8. (d) First three NOs of the pair OPDM that are not just 
localized states. The inset shows the eigenvalues of the pair OPDM. 

mass momenta that coincide with those observed in Fig. [2c) 
in the full n£ alr , evidencing the FFLO phase. 

Our results in the dilute limit are in general agreement with 
Refs.[M3 However, we have not found evidence for a fully 
paired state in the wings. According to the analysis of Ref.Uol 
[171 this phase should be stable if the magnetic field is smaller 
than half the binding energy of the bosonic pairs. This should 
occur below a critical polarization of P = 2S Z /N = 0.2 in 
the strong coupling limit. In our calculations, for all values 
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of density and magnetization considered (not shown here), we 
have only seen a fully polarized phase surrounding the core. 

In conclusion, we have numerically studied the attractive 
Hubbard model in a parabolic trapping potential, in different 
density and magnetization regimes. In the unpolarized limit, 
the system tends to be superconducting, but when the den- 
sity increases, an insulating Fock state of localized pairs forms 
in the center of the trap, displacing the superconducting state 
toward the boundaries. When the spin population is imbal- 
anced, two different limits are realized: At low densities we 
have found a structure with two regions, a well defined su- 
perconducting inner core, and fully polarized metallic wings 
that effectively behave as a ID gas of non-interacting spinless 
fermions. As the density increases, four phases emerge as a 
function of S z : a Fock state in the center of the trap is grad- 
ually replaced by free J. -fermions moving in a fully filled f 
background. A superconducting shell separates it from the 
fully polarized wings. At finite magnetizations, the super- 
conducting state becomes partially polarized, and can be de- 
scribed as a FFLO state with an oscillating order parameter 
and tightly bound pairs with a finite center-of-mass momen- 
tum. In the mean-field theory, the FFLO phase is regarded to 
be stable in a very narrow window around the Fermi surface, 7 
but we find that the presence of the trap helps to stabilize it by 
phase-separating the system. 
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